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ABSTRACT 

Deep narrowband (NB359) imaging with Subaru telescope by Iwata et al. have de- 
tected surprisingly strong Lyman continuum (LyC; ~ 900 A in the rest-frame) from 
some Lyman a emitters (LAEs) at z = 3.1. However, a possibility of redshift misiden- 
tification was not rejected with previous spectroscopy due to a narrow wavelength 
coverage. We here present results of new deep spectroscopy covering the observed 
4,000-7,000 A with VLT/VIMOS and Subaru/FOCAS of 8 LAEs detected in NB359. 
All the eight objects have only one detectable emission line around 4,970 A which is 
most likely to be Lya at z = 3.1, and thus, the objects are certainly LAEs at the 
redshift. However, five of them show a ~ 0."8 spatial offset between the Lya emission 
and the source detected in NB359. No indications of the redshifts of the NB359 sources 
are found although it is statistically difficult that all the five LAEs have a foreground 
object accounting for the NB359 flux. The rest three LAEs show no significant offset 
from the NB359 position. Therefore, we conclude that they are truly LyC emitting 
LAEs at z = 3.1. We also examine the stellar population which simultaneously ac- 
counts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet 
of the LAEs. We consider the latest statistics of Lyman limit systems to estimate the 
LyC optical depth in the IGM and an additional contribution of the bound-free LyC 
from photo-ionized nebulae to the LyC emissivity. As a result, we find that stellar 
populations with metallicity Z ^ 1/50Z© can explain the observed LyC strength only 
with a very top-heavy initial mass function (IMF; (m) ~ 50M©). However, the critical 
metallicity for such an IMF is expected to be much lower. A very young (~ 1 Myr) and 
massive (~ 100 M Q ) extremely metal-poor (Z^5x lO -4 ^) or metal-free (so-called 
Population III) stellar population can reproduce the observed LyC strength. The re- 
quired mass fraction of such 'primordial' stellar population is ~ 1-10% in total stellar 
mass of the LAEs. We also present a possible evolutionary scenario of galaxies emitting 
strong LyC and implications of the primordial stars at z ~ 3 for the metal enrichment 
in the intergalactic medium and for the ionizing background and reionization. 

Key words: cosmology: observations — galaxies: evolution — galaxies: high-redshift 
— intergalactic medium 



* This work is based on data collected with the Subaru Telescope operated by the National Astronomical Observatory of Japan 
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1 INTRODUCTION 

The first generation of stars in the universe is the stellar 
population without any metal el ements, or so-called Pop - 
ulation III (Pop III) stars (e.g., iBromm fc Larsod 120041 ). 
This pop ulation is expected to b e as massive as ~ 100 
Mq (e.g., IBromm fc Larso"nl [20041 and is thought to play 
an important role on t he cosmic reionization at z > 6 (e.g., 
lLoeb fc Barkanal 120011 ). How long did such metal- free star 
formation last in the universe? Although the metal enrich- 
ment in the intergalactic medium (IGM) is uncertain, if 
it is inefficient, metal-free hal oes may fo rm until z = 5 
(iTrenti. Stiavelli, fc Shulj 120091) , z = 2.5 (|Tornatore et all 
120071 ). or even z < 2 (| Johnson! [2O10l ). The metal mixing 
process in the interstellar medium (ISM) of a galaxy is also 
uncertain, but a significant amount of metal-fre e gas may co- 
exist with enriched gas for a few hundred Myr l|Pan fc Scald 
120071 ). If it is true, metal-free stars may continue to form in 
galaxies for such a long time. 

The He n A1640 emission line is proposed as a signature 
of Pop III stars because of their very hard ion iz ing spectrum 
expec ted theoretically (|Schaererl 120021 . 120031 ). iNagao et al.l 
(2008) made a wide and deep survey of the emission line at 
z ~ 4 and obtained an upper limit on Pop III star formation 
rate density at the redshift. On the other hand, the emission 
line has been found in a com posite spectrum of Lyman break 
galax i es (LBGs) at z ~ 3 JShaplev et al l 120031 ; iNoll et all 
l2004l ).[j imenez fc Haimanl |2006l ) claimed that it was a sig- 
nature of the metal-free stellar population in the LBGs. 
However, this He n line is broad (FWHM ~ 1500 km s" 1 ), 
and thus, it is probably caused by stell ar winds from Wolf- 
Rayet star s wit h 'normal' metallicity (IShaplev et al.l 20031 ; 



INoll et al.ll2004l ; iBrinchmann. Pettini. fc Charlotll2008l ) 



Another signature of Pop III star s is a very large Lya 
equiv alent width (EW) as > 240 A (Ma lhotra fc Rhoadsi 
120021 ). Some surveys of Lyman a emitters (L AEs) found 
galax ies with suc h a large Lya EW at z = 3.1 (Na kamural 
l20ld). z = 4.5 (iMalhotra fc Rhoadsi 120021 ) and z = 5.7 
( Shimasaku et al.l 2006t ). However, a clumpy dusty medium 



may boost the EW only apparently (|Neufeldlll99ll ). In addi- 
tion, there is a large uncertainty on the EW measurements 
because LAEs are so faint that we may not measure their 
continuum level accurately enough. 

The stellar population with metal mass fraction (or 
metallicity) Z < 10 -5 , which is < 1/2000 Z e , is 
classified as extremely metal-poor (EMP) stars (e.g., 
iBeers fc~C hristlicb 2005) and may be the second generation 
of stars l|Umeda fc Nomotd 120031 ). Hundreds of low-mass 
EMP stars have been fou nd in the halo of the Galaxy (e.g., 
IBeers fc Christlied 120051 ). Because of their long lifetime, 
these stars are survivors of the early stage of the formation 
of the Galaxy. On the other hand, it is expected that their 
high-mass counterpart exis ted in the early days and died 
out u ntil the current epoch l|Tumli nson 2006|; iKomiva et al.l 
120071 ). Yet, there is no direct observational evidence of such 
massive EMP stars at high-z. 

Although there is no confident observational signature 
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of primordial massive stellar populations such as Pop III 
and massive EMP stars so far, these populations should ex- 
ist in the early universe. These stars emit strong ionizing 
radiation which should ha ve played an important role in the 
cosmic reionization (e.g.. lLoeb fc Barkand [20011 and have 
affected the galaxy form ation in the subsequent epoch (e.g., 
ISusa fc Umemurall2004 ). Therefore, looking for these stellar 



populations is highly important to understand the feedback 
process by the first generation of stars and galaxies as well 
as to prove their existence. 

Very recently, some candidates of galaxies at z > 7 have 
been discovered by the standard drop-out technique from the 



Space Telescope ( 


Bouwcns et al. 


l2010d: lOesch et alj l20ld: 


iBunker et al.ll2009 


: McLure et al 


2010:lYan et al.ll20ld). In- 


terestingly, Bouwens et al. (2010b) found that the ultravio- 



let (UV) colour of the z > 7 candidate galaxies was so blue 
that the galaxies may be composed of EMP or metal-free 
stellar populations. However, we should confirm that the 
redshifts of the galaxies are really z > 7 with spectroscopy, 
which is too difficult to do with current facilities. 

Ilnouel l|2010l ) (hereafter Paper I) has proposed a new 
method to find the hypothetical stellar populations like mas- 
sive EMP and metal-free stars at z < 4. When Lyman con- 
tinuum (LyC) emitted by stars escapes from galaxies, the 
LyC emitted by photo-ionized nebulae around ionizing stars 
may also escape. This bound-free nebular LyC has a peak 
just below the Lyman limit, so that a spectral 'bump' ap- 
pears at the Lyman limit. We call this 'Lyman bump' or 
more precisely 'Lyman limit bump'. The strength of the Ly- 
man bump depends on the hardness of the stellar LyC. Thus, 
more metal-poor galaxies will show stronger Lyman bump. 
On the other hand, neutral hydrogen remains in the IGM 
even after the completion of the reionization and it signif- 
icantly absorbs the Ly C. It makes difficult to find LyC or 
Lyman bump at z > 5 (|lnoue fc Iwa ta 2008). 

Are there galaxies with a Lyman bump in the real 
uni verse? This paper intends to show that the answer is 
yes. Ilwata et all (|2009l ) (hereafter 109) discovered 10 LAEs 
and 7 LBGs at z ~ 3 with a significant leakage of their 
LyC captured in a deep narrowband image taken with the 
Subaru/Suprime-Cam (S-Cam). Some of the LAEs are in- 
deed brighter in LyC than in non-ionizing UV, strongly sug- 
gesting the presence of the Lyman limit bump. However, 
there was a possibility that these objects are at lower-z 
because of a narrow wavelength coverage of previous spec- 
troscopy which only confirmed an emission line. 

The rest of this paper is organized as follows. In section 
2, we first describe the sample of the possible LyC emit- 
ting LAEs discovered by 109. In section 3, we present re- 
sults of follow-up spectroscopy with VLT/VIMOS and Sub- 
aru/FOCAS. In section 4, we thoroughly examine the red- 
shift of the LAEs (i.e. the reality of the detected LyC) with 
the follow-up deep spectra and conclude that at least three 
LAEs are real LyC emitters. In section 5, we compare the 
observed strength of the LyC with the Lyman bump model 
proposed by Paper I and show that the LAEs are likely to 
contain massive metal-free or EMP stars. In section 6, we 
present a summary of our results, an evolutionary scenario 
of the LAEs, and a few implications of the existence of pri- 
mordial stars at z ~ 3. 

We adopt the AB magnitude system (lOkdl 19741) to de- 
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scribe object magnitudes and colours in section 5. The stan- 
dard flat ACDM cosmology with h = 0.7, Qm = 0.3, and 
Qa — 0.7 is adopted if it is required. 



2 THE SAMPLE OF LAES EMITTING 
STRONG LYMAN CONTINUUM 

We mainly deal with the sample of LAEs detected in their 
LyC by 109. Since the detected LyC is too strong to be 
explained with a standard stellar population model, 109 re- 
ported the detections but left their nature as a mystery. 

We started from a sample of LAEs and Lycv blobs 
(LABs) selected through a Subar u/S-Cam narrowband 
NB497 imaging in th e SSA22 field l|Havashino et alj|2004 
iMatsuda et al.l l2004 hereafter M04). The selection crite- 
ria of the LAEs were NB497 = 20.0-26.2 AB (2."0 <t>, 
> 5u) and the observed equivalent width of the emission 
line > 80 A. Among them, we had 125 galaxies with spectro- 
scopic redshift z spcc > 3.0 measured with Sub aru/FOCAS 
and Keck/DEIMOS l|Matsuda et al.ll2005l . [ioOrj . in prepara- 
tion; Yamada et al. in preparation). 109 performed a very 
deep imaging with another narrowband filter NB359 with 
Subaru/S-Cam (program-ID: S07B-010) in the same field. 
The NB359 filter exactly captures LyC (~ 880 A) for z = 3.1 
galaxies. In the NB359 image, 109 found 10 objects identi- 
fied with the spectroscopically confirmed LAEs. The iden- 
tification of NB359 (and R) sources with the spectroscopic 
LAEs was done by the following procedure: (1) object de- 
tections in NB359 and in R, (2) identification of each LAE 
with the R object nearest from the barycentre of NB497 
intensity, and (3) either identification of a NB359 object 
within l."4 from the barycentre of R intensity or > 3<r de- 
tection at NB359 aperture photometry with l."4 diameter 
around the barycentre of R intensity. Here we present re- 
sults of follow-up deep spectroscopy with VLT/VIMOS and 
Subaru/FOCAS for 5 out of the 10 LAEs in 109: the objects 
a-e in Table 1. 

In addition, there are other ~ 30 objects detected in the 
NB359 image and identified with the LAE candidates (i.e. 
not confirmed spectroscopically yet) by the same procedure 
described above. Among them, we also present results of 
deep spectroscopy with VIMOS for 4 objects: f-i in Table lQ 

In Table 1, we summarise some properties of the sam- 
ple LAEs: coordinates in NB359 (LyC), offsets and position 
angles (PAs) of R and NB497 (Lya) positions against the 
positions in NB359, and magnitudes in NB359 and R. As 
found from the table, there are small offsets in the three 
bands; the offsets between R and NB359 are less than 0."4 
and those between NB359 and NB497 are 0."1-1."0. The 
offsets of the continuum emission traced by R and NB359 of 
the LAE sample are small and not significant relative to the 
positional uncertainty between the two bands (~ 0."25; 109). 
We should note that the J?-NB359 offsets of the LBG sample 
reported in 109 is much larger and significant (~ 0."97 on 
average). On the other hand, in some cases, the offsets be- 
tween NB359 and NB497 of the LAE sample are significant. 
This point will be discussed in detail in §3.2 and §4.2. 



1 We recently found that the spectrum of the object g were taken 
with Subaru/FOCAS previously (Yamada et al. in preparation). 



3 SPECTROSCOPY OF THE LAES 

3.1 Previous medium- resolution spectroscopy 

We have several sets of previous spectroscopy with Sub- 
aru/F OCAS and Keck/DE I MOS taken in 2003, 2004, and 
2005 (|Matsuda et alj 120051 . 120061 . in preparation; Yamada 
et al. in preparation) and have confirmed that all the 10 
LAEs reported by 109 have a prominent emission line at 
around 4,970 A which was identified as Lya at z ~ 3.1. 
DEI MOS observations in 2 004 and FOCAS observations in 
2005 (Ma tsuda et aLlfeoOrj . in preparation; Yamada et al. in 
preparation) were made with a relatively high spectral res- 
olution (R ~ 2, 000) and rejected the possibility that the 
emission line was the [O n] A3727 doublet at z = 0.33. 
However, the FOCAS observations in 2005 could not ex- 
clude the possibility that the line was C IV A1549, He II 
A1640, C ill] A1909 or Mg II A2798 from an AGN at red- 
shift 2 = 2.2, 2.0, 1.6 or 0.78, respectively, or H/3 from a 
very faint galaxy at z — 0.022 because of a narrow wave- 
length coverage (AA -200 A). Since the full width at half- 
maximum (FWHM) of the line was as narrow as 200-600 
km s -1 , the objects would be a type 2 if they were AGNs. 
However, some of the objects are spatially extended in the 
NB497 — BV image (i.e. line image), making the AGN in- 
terpretation unlikely. For example, the object a is identified 
as LAB by M04 and the object c has a FWHM of the line 
intensity profile of l."3 against a l."0 PSF (a significance of 

3-CT). 



3.2 VLT/VIMOS low-resolution spectroscopy 

We performed a wide wavelength coverage (but low resolu- 
tion) follow-up spectroscopy for 5 (objects a-e) out of the 
10 LAEs reported by 109 with VLT/VIMOS (program-ID: 
081.A-0081). In addition, we took spectra of 4 (object f-i) 
LAE candidates detected in the NB359 image. Unlike the 
previous spectroscopy, we put slitlets on the positions de- 
tected in the NB359 image in order to take spectra of NB359 
emitting sources. 

The observations were made with LR-Blue/OS-Blue 
grism (R — A/ A A ~ 180 and dispersion 5.36 A/pix) and 
l."0 width slitlets on several dark nights during July to Oc- 
tober 2008. The wavelength coverage of this setting is 3,700- 
6,800 A. The plate scale is 0."205/pix. We took two masks 
(each has 4 quadrants) and the exposure time was 4 hours for 
each mask. There was an overlap area of the two masks, so 
that the object a was taken twice. The observing conditions 
were good with typical seeing of 0."8. The data reduction 
was done by a standard procedure with iraf. The flux cal- 
ibration was done with standard stars EG 274, G 158-100, 
and LTT 1020. Typical sky rms per pixel at around 5,000 A 
in the resultant two dimensional spectra is about 0.02 fj,Jy. 

Figure 1 shows close-up images around 5,000 A of the 
final two-dimensional spectra of 9 objects observed with 
VLT/VIMOS. In each panel, the area between the two 
dashed lines indicates 5 pixels (= l."0) centring around the 
position of the NB359 source. To find the NB359 source po- 
sition, we used an iraf task, GEOMATCH, between the S-Cam 
image and the VIMOS pre-image. A typical uncertainty of 
this procedure is found to be ~ 0.3 pix (= 0."06). We also 
found a few pix systematic offset between the pre-image and 
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Table 1. Coordinates of the sample galaxies. 

NB359 position Offset of R Offset of NB497 NB359 a R a Spectroscopy 

Object a (J2000) <5 (J2000) A (") PA (°) A (") PA (°) [AB] [AB] Previous Current Remarks 



a 


22:17:24.76 


+00:17:16.7 


0.23 


22 


0.38 


42 


25.84 


25.51 


FOCAS in 2005 


VIMOS/FOCAS 


109, LAB35 in M04 


b 


22:17:38.97 


+00:17:25.0 


0.26 


79 


0.85 


43 


26.28 


26.66 


FOCAS in 2005 


VIMOS/FOCAS 


109 


c 


22:17:45.87 


+00:23:19.1 


0.24 


272 


0.46 


190 


26.25 


26.72 


FOCAS in 2005 


VIMOS 


109 


d 


22:17:26.18 


+00:13:18.4 


0.21 


297 


0.88 


7 


26.68 


26.91 


DEIMOS in 2004 


VIMOS 


109 


c 


22:17:16.64 


+00:23:07.4 


0.26 


76 


1.06 


32 


26.73 


26.24 


FOCAS in 2005 


VIMOS 


109 


f 


22:17:08.04 


+00:19:31.7 


0.25 


341 


0.34 


7 


26.88 


26.52 




VIMOS 




g 


22:17:53.22 


+00:12:37.4 


0.00 




0.20 


257 


26.84 


26.84 


FOCAS in 2005 


VIMOS 




h 


22:17:29.41 


+00:06:29.0 


0.41 


144 


0.81 


189 


26.33 


26.02 




VIMOS 




i 


22:17:12.74 


+00:28:55.4 


0.16 


4 


0.12 


320 


26.26 


24.59 




VIMOS 


broad-line AGN 



a Magnitudes within a l."2 diameter aperture at NB359 or R position. 




Figure 1. Close-up around 5,000 A of the final two-dimensional spectra taken with VLT/VIMOS. The dotted lines in each panel show 
5 pixels (l."0) centring around the NB359 source position along the slitlet. North is up. 



the spectral images and corrected it based on the position of 
some bright objects detected in both images (not the LAEs 
discussed in this paper). 

As found in Figure 1, a significant spatial offset be- 
tween the emission line and the NB359 position is evident in 
some objects. The measured spatial offsets are summarised 
in Table 2 and consistent with those in Table 1 measured 
in the two narrowband images. Relatively large offset in the 
narrowband images of the object a is probably due to the 
extended Lycv emission classified as LAB. We categorise ob- 
jects a, f and g into the sample without line offset (i.e. offset 



< 0."2 = 1 pix) and objects b-e and h into the sample with 
line offset based on the spectral image. The reason why we 
rely on the spectral image is that we can measure the off- 
set between the continuum position and the pure emission 
line, whereas the NB497 intensity traces both the line and 
the continuum. The last object i is found to be a broad-line 
AGN from the spectrum. 

Figures 2 and 3 show the one-dimensional spectra ex- 
tracted at the NB359 source positions of the samples with- 
out and with the line offset, respectively. Figure 4 shows 
the spectrum of the object i which is identified as a broad- 
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Table 2. A summary of the VLT/VIMOS spectroscopy. 



Object 


Line offset (") a 


ya 


Decrement'' 


Without line offset 






a (con/emi) 


0.12 


3.100 


0.90 ±0.04 


f (con/emi) 


0.06 


3.075 


0.77 ±0.11 


cr ( fAii / pmi 1 


0.10 


3.094 


0.70 ±0.16 


composite c 




— 


0.80 ±0.07 


With line offset 








b (con) 


0.62 


— 


1.18 ± 0.29 


(cmi) 




3.090 


0.54 ±0.29 


c (con) 


0.51 




n Q4. + 015 


(emi) 




3.095 


1.38 ± 0.31 


d (con) 


1.29 




1.49 ± 0.68 


(cmi) 




3.100 


0.77 ±0.35 


e (con) 


1.00 




1.31 ± 0.44 


(emi) 




3.065 


0.69 ± 0.24 


h (con) 


0.76 




0.73 ± 0.11 


(cmi) 




3.080 


0.87 ±0.23 


composite c (con) 






1.12 ± 0.15 


composite" (emi) 






0.86 ± 0.13 


AGN 








i (con/emi) 


0.25 


3.11 d 


0.72 ± 0.03 



a Spatial offsets between the position detected in NB359 and the 
emission line. 

b Flux density ratios between 1050-1150 A and 1250-1350 A in 
the rest-frame of the Lya redshift. Uncertainties are estimated 
from the quadratic combination of mean error of the continuum 
level and sky rms noise. The Galactic dust extinction has been 
corrected. 

c Average composite based on zj^ ya . 

d This redshift was determined from metal lines. 



line AGN. In Table 2, we also list the redshifts measured 
assuming the emission line to be Lya and the flux density- 
decrements of the continua blueward and redward from the 
line. We measured the decrements both at the NB359 source 
position and at the emission line position for the sample with 
line offset. We have corrected the decrements for the dust 
ex tinction by the Milky Way based on the extinction map 
bv lSchleeel et Zl (Til). Figure 5 shows average composite 
spectra of both samples. 



3.3 Subaru/FOCAS low-resolution spectroscopy 

We have spectra of two objects a and b taken with Sub- 
aru/FOCAS (program-ID: S08B-046). The observation was 
done with the 300B grism (dispersion: 1.34 A/pix) and the 
slit width of 0."8 during the dark night on 22 September 
2008. The wavelength coverage of this observation is 4,000- 
7,200 A. The resultant spectrum of the object a is very 
consistent with those in Figure 2 but has a higher spectral 
resolution. On the other hand, we could not find any emis- 
sion line in the spectrum of the object b. This is because 
the position angle of the FOCAS slitlet (117°) was far from 
that of the emission line from the NB359 position (43° ; see 
Table 1). The position angle of the VIMOS slitlet was 0° 
(north-south) and we did capture a part of the emission line 
as shown in Figure 1. 



1.0 

0.8 
0.6 

OS 



i i i I I i i i 



I I I 





4000 4500 5000 5500 6000 6500 

Observed wavelength [A] 

Figure 2. VLT/VIMOS one-dimensional spectra of 3 objects 
without offset between the NB359 position and the emission line. 
The dotted curves are 1-cr sky noise spectra. The vertical dashed 
lines indicate the wavelengths of the emission lines detected in 
two-dimensional spectra shown in Figure 1. In each panel, we 
note the redshift if the emission line is Lya. The two thick hori- 
zontal lines with numbers in each panel show average continuum 
levels within 1050-1150 A and 1250-1350 A in the rest-frame for 
the Lya redshift. The values are not corrected for the Galactic 
extinction. The vertical shaded regions show the ranges affected 
by Earth's atmospheric lines. 
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Figure 3. Same as Figure 2 but for 5 objects with offset be- 
tween the NB359 position and the emission line. The spectra are 
extracted at the NB359 positions. 
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Figure 4. Same as Figure 2 but for the object identified as a 
broad-line AGN. The redshift was determined by metal lines in- 
dicated by vertical dashed lines. The average continuum level of 
the longer wavelength is measured within 1300—1350 A to avoid 
an effect of the broad Lya line. 
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Figure 5. Average composite spectra of (a) 3 objects without 
offset between the NB359 position and the emission line (objects 
a, f, and g), (b) 5 objects with offset (objects b, c, d, e, and 
h), and (c) the same 5 objects but at the emission line position. 
The dotted curves are l-cr uncertainty estimated from sky noise 
spectra. We averaged the one-dimensional spectra in the rest- 
frame for the redshift assuming the emission line to be Lya. The 
thick horizontal lines with numbers indicate mean flux densities 
within 1050-1150 A and 1250-1350 A. At the latter wavelength 
range, we normalised the spectra before making an average. The 
values are not corrected for the Galactic extinction. The vertical 
dashed lines indicate wavelengths of emission/absorption lines. 
The vertical shaded regions show the ranges affected by Earth's 
atmospheric lines. 



4 REALITY OF THE LYMAN CONTINUUM 

In this section, we discuss whether NB359 fluxes detected 
from the sample LAEs are truly LyC or not. Namely, we 
discuss whether the redshifts of the continuum sources are 
2 ~ 3.1 or not. As found in Figure 1, some objects show 
a spatial offset of the emission line at around 4,970 A. We 



deal with the two samples with or without the line offset 
separately. 



4.1 Sample without line offset 

4-1.1 Line identification 

From the objects a, f, and g, we clearly detect an emission 
line at around 4,970 A and no significant spatial offset from 
the position detected in NB359. Therefore, the emission line 
is highly likely to come from the continuum source detected 
in NB359 (see §4.1.2 about the probability of foreground 
contamination). Finally, we conclude that the emission line 
is Lya at z ~ 3.1 and the NB359 flux is truly LyC, based 
on the following considerations. 

We could not detect any significant emission/absorption 
lines except for the 4,970 A line from the 3 galaxies shown 
in Figure 2. Upper limits (3-cr) on the line flux ratio relative 
to the 4,970 A line are < 0.10-0.18 at 4,000 A, < 0.10-0.18 
at 4,500 A, < 0.11-0.20 at 5,000 A, < 0.15-0.24 at 5,500 A, 
< 0.20-0.35 at 6,000 A, and < 0.24-0.42 at 6,500 A. These 
limits are estimated from a comparison between a typical sky 
rms noise for 5-pix in spatial direction at a wavelength and 
the 4,970 A line flux. We have assumed that the postulated 
emission line has the same width as the 4,970 A line. This 
assumption is justified if the lines are narrower than about 
1,600 km s _1 (i.e. R ~ 180). These 3-cr upper limits of the 
line flux ratio indicate that we can exclude the existence of 
other emission lines even if the lines are a factor of 0.1 to 
0.5 weaker than the 4,970 A emission line. 

Let us identify the 4,970 A line. We consider 7 possible 
identifications of the line. If the 4,970 A line was C IV and 
the galaxies were AGNs at z ~ 2.2, we should have Lya 
around 3,900 A. Ba sed on the QSO average spectrum by 
iFrancis et alj jl991), the flux ratio of Lya to C IV is 1.6. 
Thus, we should detect the Lya in our spectra with a high 
significance. However, we could not. We could not detect He 
II and C ill] which have flux ratios rel a tive to C IV of . 3 and 
0.5, respecti vely (|Francis et al.l Il99ll ; iHeckman et al. 1995; 
lKroriklll999T ), either. If the 4,970 A line was He II or C ill] 
and the galaxies were AGNs at z ~ 2.0 or 1.6, we could de- 
tect C IV which is stronger than the two lines l|Francis et al.l 
ll99ll ; lHeckman et al.|[l995l ; iKrolikll 19991 ) at around 4,690 A 
or 4,030 A. In the Mg II case, we could detect a stronger [O 
ill emission li ne from narrow-line regions at around 6,630 A 
(|Krolikl [l999l V In the [O n] or the H/3 cases, we could de- 
tect stronger [O ill] emission lines at around 6595/6659 A 
or 5,070/5,120 A, respectively. However, we could not at all. 
The [O n] interpretation had been already ruled out by the 
previous medium-resolution spectra for the objects a-e and 
g. In the case of Lya from star- forming galaxies at z = 3.1, 
we do not expect any other stronger emission lines within 
the wavelength coverage. Therefore, the single emission line 
in the wide wavelength coverage strongly suggests that the 
line is Lya and the objects lie at z ~ 3.1. 

If the galaxies are at z = 3.1, we expect Lya decrement 
to be 0.75 ± 0.0 3 on average at around 1,100 A in rest-frame 
of th e galaxies (|lnoue fc I wata 2008; Fauche 7-Giguere et al.l 
2008) (Note that we are looking at Lya forest at z = 2.7 in 
this wavelength.). Let us see the decrement in the spectra 
of the 3 objects (a, f and g) and the AGN (object i). Com- 
paring the flux density within the rest-frame 1,050-1,150 
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A to that within 1,250-1,350 A with the assumptions of the 
4,970 A line to be Lya and a flat intrinsic spectrum between 
1,100 A and 1,300 A in the rest-frame, we find decrements of 
0.7-0.9 after correcting for the Galactic extinction as sum- 
marised in Table 2. Although the decrement of the object a 
is somewhat smaller than the expected, other three decre- 
ments are very consistent with the Lya forest measurement, 
again supporting our interpretation that they are z ~ 3.1 
objects. 

In summary, the 3 objects without line offset have a 
strong and narrow emission line around 4,970 A but do not 
have any other detectable emission/absorption lines. Their 
continuum shows a small break below the emission line and 
the amount of the decrement is consistent with that ex- 
pected from Lya forest if they are at z ~ 3.1. Therefore, 
we conclude that the emission line is redshifted Lya and the 
objects are truly LAEs at z ~ 3.1. 



4-1.2 Probability of foreground contamination 

Even when the emission line around 4,970 A is the z ~ 3.1 
Lya, the LAEs may have a foreground object on the line 
of sight towards them and the flux detected in NB359 
may come from the low -z interloper. This po ssibility is 
extensively discu ssed in IVanzella et al.l (|2010a ) (see also 
ISiana et all 120071 and 109). Let us estimate the probabil- 
ity to have a contamination in front of the 3 LAEs dis cussed 
in this subsection according to lVanzella et all | |2010ah . 

The 3 LAEs without line offset (a, f and g) are detected 
in NB359 with a magnitude of 26-27 AB (see Table 1) and 
the spatial offset between the Lya emission line and the 
continuum is less than 0."2 ( = 1 pix) (see Table 2 ). Based on 
the [/-band number count bv lNonino et all ((2009) of 130,200 
deg -2 in this magnitude range, we expect only 0.126% of 
chance coincidence of an object within 0."2 around another 
object. 

The object a is found from 10 objects reported by 109 
among 125 spectroscopic LAE sample. The 10 objects may 
include foreground contamination, but the probability that 
one object picked up from the 10 objects has a foreground 
within 0."2 is only 1.69^. The objects f and g are found 
from ~ 30 objects detected in NB359 among ~ 800 LAE 
candidates. If we take two objects from the ~ 30 sample, 
the probability that both two have a foreground within 0."2 
is only 0.1% and the probability that one has a foreground is 
still 6.5%. Therefore, the three objects a, f and g are unlikely 
to be contaminated by foreground objects. 



4-1-3 Spectroscopic properties 

We here note a few more measurements from the spec- 
tra. Table 3 is a summary of the derived properties. The 



2 If each object has a foreground with a probability p, the proba- 
bility to have k foregrounds among n objects is ( ™ ) p k (l — p) n ~ k 
iVanzella et alj l2010al l . If we make a subsample of fc) ob- 

jects including the k foregrounds and we randomly select i objects 
from the m objects, the probability to have i, k) foregrounds 
among the i objects is ( ( ) / ( ™) ■ Since k can vary from 
j to m — i + j , the total probability to have j foregrounds among i 
objects becomes ft) ^ "P)""* (?) ( T-t ) / (?) ■ 



Table 3. Spectroscopic properties of the LAEs without line offset. 



Object 




FWHM" (km s" 1 ) 


EWo' (A) 




a 


3.100 


1190 


20.4 ± 1.1 


< 0.24 


f 


3.075 


1070 


29.3 ±3.4 


< 0.40 


g 


3.094 


1130 


78.5 ± 12.3 


< 0.25 


composite^ 




1140 


43.2 ±2.9 


< 0.21 



a Full width at half-maximum of Lya. 

b Rest-frame equivalent widths of Lya. Uncertainties include only 
statistical errors in estimating the line flux based on sky rms noise. 
c 3-<t upper limits on the He II line flux relative to Lya estimated 
from sky rms noise. The line width is assumed to be the same as 
that of Lya. 
d Average composite. 



FWHMs of Lya are smaller than the velocity resolution of 
our observations (~ 1,600 km s^ 1 ). Thus, the line width 
is unresolved. The EWs of Lya in the rest-frame are 20- 
80 A in individual spectra and 43 A in the average spec- 
trum, which are not very large as expected for Pop III stars 
(e.g., iMalhotra fc Rhoadsl I200I ; ISchaereJliool ). If LyC es- 
cape fraction is relatively large, however, the expected EW 
of Lya decreases. The EW is further reduced by neutral 
hydrogen in the ISM and in the IGM because of resonant 
scattering. Therefore, the measured Lya EWs do not reject 
the existence of Pop III stars in the three LAEs. 

The He 11 A1640 line is not detected and its 3-tr up- 
per limits relative to Lya flux (-Fkeii/^Lyar) are < 0.2-0.4, 
where we have assumed tha t the He 1 1 line has the same 
width as Lya. According to Schaerer (200j|), the intrinsic 
fkeii/ 'Fi,ya < 0.02 is expected for a Pop III star cluster. 
Since only Lya photons will be scattered resonantly in the 
surrounding ISM and IGM, the observed FHell/-FLy Q can be 
enhanced significantly. On the other hand, the He 11 line is 
detectable only in a short time-scale (< a few Myr) from 
the Pop III star formation ((Sc haerer 20021). In any case, the 
obtained upper limits on Fn e ii/Fi,y a are not strong enough 
to reject the presence of Pop III stars. 

We could not detect any other emission/absorption lines 
in the spectra. The 3-ct flux upper limits of emission lines 
are 0.1-0.4 relative to Lya depending on the wavelength 
(§4.1.1). In this respect, the three LAEs presented here are 
different from the Lynx arc at z — 3.357 and a peculiar LAE 
at z — 5.563 emitting strong metal emission lines such as N 
ivl A1486, C iv A1549, O ml A1666 and C ml A1909 as well as 



Lya ( Fosburv et al.l [20031; [Raiter. Fosburv. fc Teimoorinial 
|2010| ; iVanzella et allhoiObl ). The observed fluxes of N iv] 



and C iv, which are in the wavelength coverage of our spec- 
tra, relative to Lya are 0.1-0.3, which could be marginally 
detected in our spectra if there were. These two peculiar 
galaxies are likely to have very hot (~ 10 5 K) exciting stars 
i|Fosburv et al.ll2003l; (Raiter. Fosburv. fc Teimoorinial l2oioh 
or small AGN ( Vanzella et al.l 2010bl ) and are as low metal- 
licity as Z = 1/20^0. On the other hand, our three LAEs 
may have massive EMP or Pop III stars, whose effective 
temperatures are also ~ 10 5 K, as shown later (§5). Never- 
theless, we could not detect such metal emission lines. This 
may imply that the metallicity of our three LAEs is much 
lower than Z — 1/20Zq which is consistent with the pres- 
ence of EMP or Pop III stars. 
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4.2 Sample with line offset 

4-2.1 Line identification 

We detected a prominent emission line at 4,970 A around 
the objects b-e, and h as shown in Figure 1. However, the 
line offsets from the NB359 source position by 0."5-l."3, 
which correspond to 3.8-9.9 kpc (proper) if the objects lie 
at z = 3.1. In the emission line positions, we find no other 
emission/absorption lines. After a very similar discussion to 
§4.1.1, we conclude that the emission line is likely to be Lya 
at z ~ 3.1. In the NB359 positions, we do not find any 
significant emission/absorption lines, either. 

We measured Lya decrements in both NB359 and emis- 
sion line positions of each object, assuming the emission line 
to be Lya, and the results are summarised in Table 2. We 
detect a significant (> 2-<r) decrement at the continuum po- 
sition of the object h and marginal (1-2-cx) ones at the line 
position of the objects b and e. Other cases are not conclu- 
sive because of not-significant continua although some cases 
may imply no decrements. The decrements in the compos- 
ite spectra are not conclusive, either, although that at the 
NB359 positions exceeds unity, indicating contaminations of 
foreground objects. 

In summary, the emission lines detected in our spectra 
are likely to be Lya, and thus, these objects are LAEs at z ~ 
3.1. However, almost no information about redshift of the 
NB359 continuum sources was extracted from our spectra. 
These sources may be physically associated with the close 
LAEs and their redshift may be z ~ 3.1. Or these are just 
foreground objects apparently close to the LAEs. 



4-2.2 Probability of foreground contamination 

As discussed in §4.1.2, there is a possibility that a foreground 
interloper accounts for the NB359 flux (and other optical 
band fluxes). The 5 objects we are discussing here show an 
offset from Lya line by 0."8 on average (see Table 2). Let us 
estimate the probability of such fo reground contamination 
according to lVanzella et all (|2010al l as done in §4.1.2. 

The observed NB359 magnitudes of the 5 objects (b- 
e and h) are 26-27 AB same as the 3 LAEs without line 
offset. Then, the probability to have an object within 0."8 
around another object is 2.02% based on the spatial densit y 
of U = 26-27 AB sources reported bv lNonino et all |200g| ). 

The objects b-e are selected from 10 LAEs which were 
found from 125 spectroscopic LAEs and reported by 109. 
The probability that all the four taken from the 10 sample 
have a foreground within 0."8 is only 0.8%, but the prob- 
abilities that 3, 2, 1, or objects have a foreground are 
5.7%, 20.6%, 39.8%, or 33.2%. The object h are found from 
~ 30 objects detected in NB359 among ~ 800 LAE candi- 
dates. The probability that the object is contaminated by 
a foreground within 0."8 is 53.8%. Therefore, a few of the 
five objects are likely to be contaminated by a foreground. 
However, it is very difficult to conclude that all the five are 
contaminated. 

As found in Table 2, the objects d and e have a rela- 
tively large offset (> l."0) and no indication of Lya decre- 
ment at the NB359 positions. It may suggest that these are 
contaminations. On the other hand, the object c has the 
smallest offset and a possible indication of decrement, may 



suggesting that it is a real LyC emitter. However, it is im- 
possible to conclude which objects are contaminations using 
the current data. 



4-2.3 Possible causes of Lya line offset 

We can still expect that a few objects among the 5 are as- 
sociated with LAEs at z ~ 3.1 and the NB359 traces the 
LyC according to the statistical argument obtained in §4.2.2. 
In other words, the offset between the Lya emission line 
and the continuum could be real in a few cases. For ex- 
ample, it is possible that a galaxy is composed of multiple 
substructures and every component emits both the contin- 
uum and the emission line. However, the line flux relative 
to the continuum (i.e. EW) are different from each other. 
In fact, s uch a situation is observed in local galaxies emit- 
ting Lya [jOstlin et alj|2009h . Then, if we observe the galaxy 
as a single object with an insufficient spatial resolution, the 
barycentre of the emission line intensity and of the contin- 
uum will offset. 

Another possibility is that Lya emission comes 
from 'cold accretion' into galaxies JPiikstra fc Loebl 12009] ; 
iGoerdt et aUl20ld ; IFaucher-Giguere et al.ll2010l ). The 'cold 
accretion (or stream) ' is thought to be the main mode to 
feed gas to galaxies |Pekel et al.1 2009), and its tempera- 
ture is ~ 10 4 K, so that it emits Lya. Galaxies have mul- 
tiple streams a nd Lya emission from the streams is spa- 
tially extended jDiikstra fc Loebl 120091 ; [Goerdt et al.ll20ld ; 
IFaucher-Giguere et al. 2010T I. If we observe such an ex- 
tended Lya emission with insufficient spatial resolution, the 
barycentre of the emission may offset from the stellar com- 
ponent. 

Therefore, the Lya offset may not be a surprising event 
for LAEs in general. However, we have no further clue to 
resolve which is a real z ~ 3.1 object for the moment. 



5 REST UV TWO-COLOUR DIAGRAM 

In this section, we try to interpret atypical UV colours of 
the LAEs (and also LBGs) detected in LyC. Figure 6 shows 
a two-colour diagram in the rest-frame UV of the LAEs and 
LBGs. We consider that the 3 LAEs without line offset (filled 
circles) are emitting LyC and a few of the 5 LAEs with line 
offset (open circles) are possibly emitting LyC. The 7 LBGs 
reported by 109 are also shown by filled squares. Addition- 
ally, we show the AGN (object i) found in this study. The 
colours are measured at the barycentre in R with a circu- 
lar aperture whose diameter is chosen from l."2 to 2. "6 for 
each object so as to include its total flux as done in 1090 
Note that there are differences from Table 1 because of the 
different method of measurements. 



3 There was a revision of zero-points in V, R, and i' after publi- 
cation of 109 and the resultant shifts in the two colour diagram as 
follows: A(NB359 — R) = +0.1 and A(V - i') = +0.1. Although 
the zero-point revision is almost settled, we conservatively apply 
additional systematic uncertainties of 0.1 mag to the error-bars of 
the observed points by the quadratic combination with statistical 
uncertainties in photometry. 
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Figure 6. Rest UV two-colour diagram for z ~ 3.1 LyC emitting galaxies. The vertical axis, V — i', indicates non-ionizing UV spectral 
slope, and the horizontal axis, NB359— R, indicates LyC-to-UV flux density ratio. The filled circles with errorbars are the 3 LAEs without 
line offset and the open circles with errorbars are the 5 LAEs with line offset. The squares with errorbars are LBGs reported in Iwata et 
al. (2009). The plus mark with errorbars is an AGN emitting LyC found in this study. The intrinsic colours of stellar population models 
are shown by diamond (the model D in Table 4), asterisk (C), x-marks (large: Bl; small: B2), triangles (large: Al; small: A2) and inverse 
triangles (large: Acl with 100 Myr; small: Ac2 with 100 Myr). The two horizontal arrows show IGM minimum and median attenuations 
for z = 3.1. The two diagonal arrows show dust attenuations with E(B — V) = 0.1 for the SMC law (solid arrow) and the Calzetti law 
(dashed arrow). The shaded regions indicate the regions explained by the stellar population model Al with a combination of IGM and 
dust attenuations. 



Table 4. Models of stellar population and SED. 



Model 


Z/Zq 


P 




"Mow 


(m) 


SF history 


age 


SED reference 


Al 


1/50 


-2.35 


IOOMq 


1M Q 


3.1M 


Instantaneous 


1 Myr 


STARBURST99 (v. 5.1) 


A2 


1/5 


-2.35 


100M q 


1M q 


3.1M 


Instantaneous 


1 Myr 


STARBURST99 (v. 5.1) 


Acl 


1/50 


-2.35 


100M Q 


1M 


3.1M 


Constant 


10, 100, 1000 Myr 


STARBURST99 (v.5.1) 


Ac2 


1/5 


-2.35 


IOOMq 


1M 


3.1M 


Constant 


10, 100, 1000 Myr 


STARBURST99 (v.5.1) 


Bl 


1/50 


-0.1 


IOOMq 


1M 


48M 


Instantaneous 


1 Myr 


STARBURST99 (v.5.1) 


B2 


1/5 


-0.1 


IOOMq 


1M 


48M 


Instantaneous 


1 Myr 


STARBURST99 (v.5.1) 


C 


1/2000 


-2.35 


500Mq 


5OM 


IIIMq 


Instantaneous 


1 Myr 


Schaerer (2003) 


D 





-2.35 


5OOM 


5OM 


IIIMq 


Instantaneous 


1 Myr 


Schaerer (2003) 



5.1 Stellar population models 

In order to interpret the UV colours of the LyC emitting 
LAEs, we construct models of spectral energy distributions 
(SEDs) of galaxies. The SED of pure stellar populations de- 
pends on metallicity, initial mass function (IMF), star for- 
mation history, and age. Since we are dealing with LAEs 
(and LBGs), we consider sub-solar metallicities: Z = 1/5 Zq 
and 1/50 Zq. In addition, we consider two more metallici- 
ties: Z = 1/2000 Z e (EMP) and (metal- free). SEDs for the 
former two metallicities are generated by the population syn - 
thesis code starburst99 version 5.1 jLeitherer et al.|[l999h . 
Thos e for the latter two are taken from ISchaeren 1 20021 . 
2003). T he IM F is basically assumed to be Salpeter's one 
l|Salpeter|[l955l ): 4>(m)dm oc mPdm with p — —2.35. Ad- 



ditionally, we consider an extremely top-heavy case with 
p — —0.1. The mass range of the IMF is assumed to be 
m up = 100 Mq and mi ow = 1 Mq for Z — 1/5 Zq and 
1/50 Z cases but to be m up = 500 Mq and mi ow = 50 
Mq for the EMP and metal-free cases. The average stel- 
lar masses are 3.1 Mq, 48 Mq, or 111 Mq depending on 
the IMF. We consider age of 1 Myr after an instantaneous 
star formation or 10 Myr to 1 Gyr constant star formation. 
Table 4 is a summary of the models and their parameters. 

In Figure 6, the intrinsic (i.e. no dust and IGM atten- 
uations) colours of the stellar population models are shown 
by diamond (metal-free: the model D), asterisk (EMP: the 
model C), x-marks (extremely top-heavy IMF but normal 
sub-solar metallicity: the models Bl and B2), triangles (nor- 
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mal IMF and sub-solar metallicity: the models Al and A2), 
and inverse triangles (normal IMF and sub-solar metallic- 
ity but constant star formation: the models Acl and Ac2). 
NB359— R colours of the 3 LAEs without line offset (large 
filled circles) and LBGs (small squares) seem to be explained 
by a model with normal IMF and metallicity (models A and 
Ac). However, the IGM attenuation makes NB359— R redder 
as described in the next subsection. 



5.2 IGM attenuation 

Since the IGM attenuation, especially for the LyC, has 
a stochastic nature, we e stimate the amount b ased on a 
Monte-Carlo simulation bv llnoue fe Iwa ta (2008) (hereafter 
1108). These authors assumed an empirical distribution func- 
tion of the intervening absorbers in redshift, column density, 
and Doppler parameter spaces which was derived from the 
latest observational statistics in that time. Their simulation 
reproduces the Lya depressions at z = 0-6 very well. 

The LyC opacity is mainly determined by the Ly- 
man limit systems (LLSs; clouds with log 10 (A r Hi/cm -2 ) ^ 
17.2), we compare the number density evolution adopted 
in 1108 with new survey results after the paper in Fig- 
ure 7. We find that o bserva tional data are not converged 
yet; iProchaska et al.l (|20 ldl ) obt ained smaller numb er of 
LLSs at z < 4 t han t hose by iPeroux et al.l (|2005l ) and 
ISongaila fc Cowid i|2010l ). We also find that the simulation 
by 1108 adopted a number density tracing upper bounds of 
observations. In order to match the new observations, we 
have reduced the LLS number density, with almost keeping 
the number densities of Lya forest and damped Lya sys- 
temfl The updated number density evolution is shown by 
diamonds in Figure 7. The new number density of LLSs at 
z = 2.9 which attenuate the NB359 flux of z = 3.09 sources 
is about 60% of the previous one. 

Figure 8 shows the cumulative probability of the IGM 
attenuation through the NB359 for sources at z = 3.09 
which is the typical redshift of our LAEs (and LBGs). To es- 
timate the attenuation, we have assumed a constant source 
spectrum in f v unit. The effect of the different shape of the 
spectrum is negligible because the wavelength coverage in 
the NB359 is narrow. As found in the figure, the updated 
version (solid line) expects much smaller attenuation than 
the original 1108 (dashed line). The median attenuation is 
reduced from 1.08 mag to 0.60 mag. Thus, the LLS num- 
ber density has a significant impact on the expected IGM 
attenuation. We also find that the attenuation of the small- 
est 0.15% (e.g., smaller 3-er excess for Gaussian) is 0.35 mag 
(1108) or 0.28 mag (updated). These attenuations are real- 
ized on a line of sight without LLS at z ~ 3 and produced 
by numerous and unavoidable Lya forest. Because the up- 
date in the number density of absorbers is small for Lya 
forest, the reduction of the smallest attenuations is small. 

4 This is realized only if we change three parameters among ten 
in the empirical distribution function of absorbers. The change 
is (A,Pi,/3 2 ) = (400,1.6,1.3) to (500,1.7,1.2), where A is the 
normalisation of the number of absorbers, and /3i and 02 are the 
indices of the double power-law for the column density distribu- 
tion. Note that we should keep the number density of Lya forest 
because the Lya decrements by 1108 excellently agree with obser- 
vations. 
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Figure 7. The number density evolution of Lyman limit systems 
which have log 1 Q(A r Hi/cm -2 ) ^ 17.2. Observational data are the 
filled circles (Peroux et al. 2005), squares (Songaila & Cowie 
2010), and triangles (Prochaska et al. 2010; a factor of 1.1 larger 
than their Table 4 which are the density for log 1 Q(A r Hi/cm -2 ) ^ 
17.5). Monte-Carlo simulations are shown by the open circles (In- 
oue & Iwata 2008) and the diamonds (updated version). 
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Figure 8. The cumulative probability to have an IGM attenua- 
tion in NB359 for z = 3.09 sources smaller than that in the hor- 
izontal axis, based on our Monte-Carlo simulations: the dashed 
line (Inoue & Iwata 2008) and the solid line (updated version). 
The number of realizations of lines of sight is 10,000. The nar- 
rowband filter NB359 captures 852-905 A in the rest-frame of the 
sources. The inset shows a close-up of the smallest attenuation. 



We adopt the attenuation at the cumulative probability of 
0.15% as the minimum IGM attenuation in the following. 

In Figure 6, we show the minimum and median IGM 
attenuations by short and long horizontal arrows. We note 
that the IGM attenuations for other bands in Figure 6 are 
zero (for R and i') or negligibly small (for V). Table 5 gives 
a summary of the IGM attenuation. 



5.3 Dust attenuation 

If galaxies contain dust, they are reddened. Let us as- 
sume two types of d ust attenuation law: the Calzetti law 
l|Calzetti et all l2OO0h and the extinction law of the Small 
Magellanic Cloud (SMC). The former attenuation law was 
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Table 5. A summary of reddenings by IGM and dust for a source 
at 2 = 3.09. 





Colour excess 


£(NB359 - R) 


E(V - i') 


IGM 


Smallest 0.15% 


0.28 




(updated) 


Median 


0.60 




Dust" 


Calzetti 


0.90 


0.20 




SMC 


0.20 


0.84 



"These colour excesses correspond to E(B — V) = 0.1 in the 
rest- frame. 



deri ved from spectra of local UV-selected starburst galax- 
ies l|Calzetti et al J 1200(f ) and is routinely assumed in stud- 
ies of high-z galaxies. However, recent studies for z — 2- 
3 LBGs suggest that a steeper attenu ation law than the 
Calzetti one is favourable for some cases (|Reddv et alJ feoOG; 
ISiana et al.ll2009T ). Thus, we also adopt the SMC law which is 
the steepest among th e known attenuation/ex t inctio n laws. 
We use the model by IWeingartner &: Draind l|200ll ) which 
reproduces the empirical SMC extinction law very well. 

The Calzetti law is oris inally only for A ^ 1200 A 
jCalzetti et al.ll2000h . However, we simply extrapolate the 
original formula for A < 1200 A. Thi s could be justified 
by the study of iLeitherer et ail (|2002h up to A = 970 A 
(but see also Buat et al. 2002 ) and al so be justified by the 
model o f lWeingartner fc Draind (|20Q]J) which does not show 
any breaks or feature s up to about 700 A for graphite and 
astronomical silicate (|Draine fc Ledfl984h . Table 5 gives a 
summary of the dust reddenings. 

In Figure 6, we show the region which can be explained 
by a model with normal sub-solar metallicity and IMF and 
with a combination of dust and IGM attenuations. While a 
half of LBGs are found in the region, all the LAEs are out 
of it. We need a new model to explain these LAEs. 



5.4 Models with escaping nebular Lyman 
continuum 

Paper I has proposed the importance of an additional contri- 
bution by escaping nebular LyC. For an escape of the stel- 
lar LyC, the ionized nebulae should be 'matter-bounded', 
at least along some lines of sight, from which we can also 
expect an escape of nebular LyC produced by the recom- 
bination process. Paper I shows that this escaping nebular 
LyC makes a peaky spectral feature just below the Lyman 
limit, which we call Lyman limit 'bump' (see Figs. 3 and 4 
in Paper I). In fact, our NB359 filter exactly captures this 
Lyman limit bump at z ~ 3.1. Let us compare this scenario 
with the observed strength of LyC. 

Figure 9 shows how the Lyman limit bump blues 
NB359— R colour; the colour becomes bluer when a frac- 
tion of LyC is absorbed by nebulae and reemitted in the 
Lyman limit bump than when all the LyC escapes. Indeed, 
the bluest colour is realized when about 40% LyC is ab- 
sorbed (i.e. the escape fraction is about 60%), irrespective 
of the stellar population model. Based on Figure 9, we find 
that moderate sub-solar metallicity and normal IMF (solid 
line: model Al) cannot reach NB359-.R < 0.4 (0.7) for the 
minimum (median) IGM attenuation even if we consider the 
additional contribution by nebular LyC. On the other hand, 
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Figure 9. NB359— R colour for z = 3.1 galaxies, which corre- 
sponds to the strength of the Lyman limit 'bump', as a function 
of the escape fraction of stellar LyC: (a) minimum IGM attenu- 
ation case and (b) median IGM attenuation case. The curves are 
the models with escaping nebular LyC: the long-dashed curve for 
the model D in Table 4, the short-dashed curve for the model C, 
the dot-dashed curve for the model Bl, and the solid curve for the 
model Al. The nebular gas temperature is assumed to be 1 X 10 4 
K. 



some 2 ~ 3.1 LAEs detected in LyC show a much bluer 
colour. Furthermore, even the LAEs with NB359-i? > 0.4 
are difficult to be explained because of their red colour in 
V — i' which expect much redder NB359— R as shown in 
Figure 6. 

Figures 10-12 show the two-colour diagrams same as 
Figure 6 but with escaping nebular LyC models. In each 
panel, the curve with symbols shows the colour sequence 
as a function of the escape fraction of the stellar LyC. We 
have assumed the minimum IGM attenuation for the curves. 
The symbols indicate the colours when the escape fraction 
is 1.0 (i.e. pure stellar colour + minimum IGM), 0.1, or 
0.01. As the escape fraction decreases, NB359— R colour first 
becomes bluer than the stellar one due to the Lyman limit 
bump as shown in Figure 9, and then, the colour turns over 
when the escape fraction is about 0.6 and becomes redder 
and redder after the bluest point. In the same time, V — i' 
colour becomes redder than the stellar one because of the 
nebular continuum by two-photon and bound-free processes. 
The amount of these colour changes depends on the strength 
of the stellar LyC: larger changes by stronger LyC. 

The shaded region in each panel of Figures 10-12 shows 
the colours explained by a combination of IGM and dust at- 
tenuations for a specific stellar population model. The thick- 
ness of the shades indicate the cumulative probability of the 
IGM attenuation: higher probability, thicker. We have as- 
sumed the SMC extinction law for the dust attenuation but 
the Calzetti law cases are included in the shaded regions like 
in Figure 6. 

Figure 10 shows the cases of the normal stellar popu- 
lation with a constant star formation (Acl in Table 4). We 
have assumed Z — 1/50 Zq, very low, but, sometimes ob- 
served metallicity (hereafter we call it normally sub-solar). 
The case with Z = 1/5 Zq is always redder in NB359— R 
than the case with Z = 1/50 Zq . We find that all the LAEs 
(filled and open circles) and half of the LBGs (squares) can- 
not be explained even with the duration of the star for- 
mation of 10-100 Myr which is younger th an a typical age 
of LBGs (~ 300 Myr; IShaplev et althoOSl) . Therefore, we 
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Figure 10. Same as Fig. 6 but a comparison with the escaping 
nebular LyC scenario of a continuous star-forming galaxy with 
Z = 1/50 Zq (model Acl): (a) age from the onset of star forma- 
tion of 100 Myr and (b) 10 Myr. The dotted curves indicate colour 
sequences as a function of the escape fraction of stellar LyC with 
the minimum IGM attenuation. The symbols indicate the po- 
sitions with the escape fraction of 1.0, 0.1, or 0.01. The shaded 
regions indicate the regions explained with a combination of IGM 
and dust attenuations. From thin to thick, the cumulative proba- 
bility to have the amount of the IGM attenuation increases. The 
SMC extinction law is adopted for the dust attenuation as shown 
by the solid arrow. The nebular gas temperature is assumed to 
be 1 X 10 4 K. 



conclude that the LyC emitting LAEs should have a stel- 
lar population much younger than 10 Myr or much more 
massive than the standard Salpeter IMF. 

Figure 11 shows the cases of very young (age of 1 Myr) 
stellar populations with Z = 1 /50 Zq , normally sub-solar 
metallicity. The case with Z = 1/5 Zq is about 0.2 mag 
redder in NB359— R than the case. Assuming the standard 
Salpeter IMF (panel [a]; model Al), we find that it is still 
difficult to explain all the LAEs. On the other hand, if we 
assume an extremely top-heavy IMF whose mean mass is 
about 50 Mq (panel [b]; model Bl), the LAEs without 
line offset (filled circles) can be explained with an atten- 
uation smaller than the median. However, it is still diffi- 
cult to explain the LAEs with line offset (open circles). We 
may have a few true LyC emitters in the 5 objects because 
we have statistically rejected the possibility that all the 5 
have a faint foreground object accounting for the NB359 
flux in §4.2.2. There is another difficulty of the model Bl. 
The critical metallicity for the IMF change from the stan- 
dard to top-heavy, Z CI < 1O~ 3 Z |Bromm fc Loebl 120031 ; 
ISchneider et al.l 120031 . 120061 ) is expected to be much lower 
than Z = 1/50 Zq which is assumed in the model Bl. Thus, 
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Figure 11. Same as Fig. 10 but the cases of a very young (1 Myr) 
stellar population with Z = 1/50 Zq: (a) normal IMF (model Al; 
solid line with triangles) and (b) extremely top-heavy IMF (model 
Bl: mean stellar mass of ~ 50 Mq; dot-dashed line with crosses). 



a very massive IMF under the relatively 'high' metallicity 
in the model Bl is unlikely. Then, we conclude that it is 
difficult to explain the strength of LyC of the LAEs with 
normally sub-solar metallicity unless the critical metallicity 
for the IMF change is much higher than expected in the 
literature. 

Figure 12 shows the cases of very young (age of 1 
Myr) and massive (~ 100 Mq) stellar populations with 
Z = 1/2000 Zq (EMP; model C; panel [a]) and Z = (Pop 
III; model D; panel [b]). In these cases, we can expect a 
very massive IMF according to the critical metallic i ty sce - 
nario (Br omm fc Loebl 120031 : ISchneider etafl I2003L 120061 ). 
With the model C, we can easily explain the three LAEs 
which are confirmed to be true LyC emitters, whereas the 
LAEs with line offset are still difficult, except for the ob- 
ject b. The Pop III case (model D) can explain the objects 
b and c, and marginally explain e and h. The object d is 
still difficult but this object shows the largest offset of the 
Lycv (Table 2), and thus, the foreground probability is the 
highest. On the other hand, the object c shows the smallest 
offset and is the most likely object as a real LyC emitter 
among the 5 offset LAEs. Enough interestingly, the Pop III 
model can explain the object with a 1-20% smallest IGM 
attenuation. In summary, very young (~ 1 Myr) and mas- 
sive (~ 100 Mq) EMP stars can reproduce the observed 
colours of the three LAEs without line offset and Pop III 
stars can reproduce even the colours of most of the LAEs 
with line offset. Because a few line-offset LAEs are likely to 
be real LyC emitters, we conclude that the Pop III model is 
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Figure 12. Same as Fig. 10 but the cases of a very young (1 
Myr) and massive (~ 100 Mg) stellar population with (a) Z = 
1/2000 Zq (model C; short-dashed line with asterisks) and (b) 
Z = (model D; long-dashed line with diamonds). 

the most favourable stellar population for the LyC emitting 
LAEs. 

5.5 Models with two stellar populations 

The previous subsection shows that the LAEs detected in 
LyC by 109 are likely to contain a primordial stellar pop- 
ulation such as massive Pop III or EMP stars. But how 
much amount of these exotic stars are required in them? 
Also, are Pop III stars really compatible with a slight dust 
attenuation which is required to explain some LAEs with 
line offset and red in V — i' in Figure 12? To discuss these 
questions, let us consider a system in which a primordial 
stellar population and a normal stellar population with sub- 
solar metallicity and Salpeter IMF coexist. We assume that 
the normal population makes stars with a constant rate; 
the model Ac2 in Table 4 is adopted for this population. 
For the primordial population, we adopt the model D (Pop 
III). The stellar mass of the normal population is defined 
as the multiple of the star formation rate and the duration; 
we neglect the loss of the stellar death for simplicity. The 
mass fraction of the primordial population in the total stel- 
lar mass is a parameter for the mixture. The contribution 
of nebular continuum in LyC and in other wavelengths is 
taken into account for both populations but with different 
escape fractions: 0.5 for the primordial population and 0.01 
for the normal population^ The normal population may ex- 

5 In addition to the photo-ionized gas, we can consider the free- 
free radiation from hot gas produced by SNe and stellar winds. 
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Figure 13. Same as Fig. 6 but comparisons with models of two 
stellar populations: the cases with models D (Pop III) and Ac2 
(underlying population with a sub-solar metallicity, 1/5 Zq, and 
a normal IMF) in Table 4. A constant star formation of 100 Myr is 
assumed for the underlying population. The SMC extinction law 
is assumed for the panel (a) but the Calzetti attenuation law is as- 
sumed for (b). The median IGM attenuation is assumed for both 
panels. The short-dashed curves are the sequences of the colour 
as a function of the mass fraction of Pop III stars for dust-free 
underlying stellar population. The solid, dot-dashed, and long- 
dashed curves are the same sequences but for dusty underlying 
population with E(B — V) = 0.1, 0.2, and 0.3, respectively. Note 
that the Pop III stars are assumed to be always dust-free. The 
positions for the mass fraction of 0.1, 0.01, 0.001, and 0.0001 are 
indicated by dotted curves with squares, triangles, asterisks, and 
plus-mark, respectively. For Pop III stars, the contribution of es- 
caping nebular LyC is taken into account, assuming the escape 
fraction of 0.5. The colour sequence with other escape fractions 
is shown by the triple-dot-dashed curve as in Fig. 12. For the 
underlying population, the escape fraction of 0.01 is assumed. 

ist in dusty environment, but the primordial population are 
likely to be dust-free. We adopt two dust attenuation laws: 
the SMC extinction law and the Calzetti attenuation law as 
described in §5.3 only for the normal underlying population. 

The required mass fraction of the primordial stellar pop- 
ulation depends on the duration of the star formation in the 
underlying normal population. This is because the mass of 
the primordial population which is required to explain blue 



We have assumed that the total energy of the hot gas radia- 
tion to be 5% and 1% of the bolometric stellar luminosities for 
the primordial p opulation (as done in Paper I) and for the nor- 
mal population (Lcit herer et al,lll999h . respectively. The hot gas 
temperature is assumed to be 1 X 10 s K. However, this hot gas 
contribution is negligible as shown in Fig. 3 of Paper I. 
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NB359— R colour does not change very much because the 
colour is determined mainly by the primordial population, 
while the mass of the underlying population is proportional 
to the duration. For example, we consider 100 Myr as a du- 
ration for z — 3 galaxies (e.g., IShaplev et al.l [2003). The 
results are shown in Figure 13. For other durations, the pri- 
mordial mass fraction is roughly estimated by an inverse 
relation to the durations of the underlying star formation. 
Although we assumed the median IGM attenuation in Fig- 
ure 13, the readers can shift all the model curves along the 
horizontal axis for other IGM attenuations. 

The LAEs without offset (objects a, f and g; filled cir- 
cles in Fig. 13) require a primordial mass fraction of 0.1-10% 
depending on the IGM and dust attenuations as well as the 
star formation duration. The dust attenuation for the un- 
derlying population is modest as E(B — V) — 0-0.1 for the 
SMC case or up to 0.3 for t he Calzetti case, wh ich is reason- 
able for z ~ 3 LAEs (e.g.. iGawiser et al.|[2006l l. The bluest 
LAEs (objects b and c; open circles) require a primordial 
mass fraction of more than 10%, but we could not reject 
the possibility that their NB359 flux was foreground con- 
tamination individually (§4.2). The LAEs which are red in 
V — i' (objects d, e and h; open circles) can be explained 
by a model with the SMC law and a small amount of atten- 
uation as E(B — V) — 0.2-0.3 and with a primordial mass 
fraction of a few 0.1%. When the Calzetti law is adopted, 
the red LAEs require a large amount of dust attenuation as 
E(B — V) > 0.5 which may be too large. Note that the pos- 
sibility which the NB359 flux of these objects, especially d 
and e, are foreground contamination is the largest because 
of the largest line offset (§4.2). 

The LBGs from which 109 reported detections of LyC 
(filled squares) are found in the range of the primordial mass 
fraction of 0.01-1 % depending on the dust attenuations. 
Thus, only a very small amount of primordial stars can ac- 
count for the observed LyC of these galaxies. On the other 
hand, some of these LBGs can be explained by a normal 
population with a large escape fraction of ~ 0.5 as shown in 
Figures 10 and 11. Therefore, the primordial stars are not 
mandatory for them. 



6 SUMMARY AND DISCUSSIONS 

We have investigated the nature of the LAEs at z ~ 3.1 
detected in our deep narrowband NB359 imaging with 
Subaru/S-Cam. The NB359 captures LyC from z > 3. Thus, 
we detect LyC from the LAEs if they are at z ~ 3.1. 
These LAEs are special because of their surprisingly strong 
LyC relative to non-ionizing UV (109). Deep follow-up spec- 
troscopy with VLT/VIMOS and Subaru/FOCAS for 8 such 
LAEs presented in §3 shows that at least three of them are 
highly likely to be at z ~ 3.1 and the NB359 captures truly 
their LyC (§4.1.1 and §4.1.2). From the spectra of the three 
LAEs, we have derived the rest-frame EW of Lya which are 
not very large as expected from primordial stellar population 
models such as Pop III stars. We have also derived upper 
limits on the He II A1640 emission line which are not strict 
enough to reject Pop III star formation (Table 3). Other five 
out of the 8 LAEs show a ~ 0."8 offset between the Lya 
emission and the continuum detected in NB359 (Figure 1). 
For these LAEs, we could not determine the redshifts of the 



continuum sources (§4.2.1). However, it is statistically diffi- 
cult that all the five have foreground contamination which 
explain the NB359 detections (§4.2.2). Thus, NB359 may 
capture LyC from a few of the five LAEs with line offset. 

Very interestingly, all the LAEs reported by 109 are too 
bright in LyC to be explained by normal stellar population 
models (Figure 6). Although we reduced IGM opacity to 
match the latest LLS statistics (Figures 7 and 8), compared 
to the previous estimate, this is true even for the three LAEs 
which are confirmed to be at z ~ 3.1. Unlike traditional 
spectral models of galaxies, Paper I proposed a new model 
taking an escape of nebular LyC into account. This model 
expects a strong flux excess just below the Lyman limit: we 
call it Lyman limit 'bump'. We tried to explain the observed 
strength in LyC of the LAEs with this new Lyman 'bump' 
model. As a result, we have found that the colour of the 
LAEs cannot be explained by the stellar population with 
the standard Salpeter IMF and normally sub-solar metallic- 
ity (Z ^ 1/50 Zq) even with the Lyman 'bump' (Figures 10 
and 11a). If we assume a very young (~ 1 Myr) stellar popu- 
lation with an extremely top-heavy IMF (mean mass of ~ 50 
Mq) and Z = 1/50 Zq, the three LAEs without line offset 
can be explained with a relatively smaller IGM attenuation 
but other five LAEs with line offset cannot be explained yet 
(Figure lib). Because (1) a few of the line-offset LAEs are 
possibly to be real LyC emitters and (2) the expected criti- 
cal metallicity for a top-heavy IMF is much lower than the 
assumed Z — 1/50 Zq, we reject this model. This suggests 
that the LAEs contain a significant amount of 'primordial' 
stars such as EMP or Pop III stars. 

A very young (~ 1 Myr) and massive (~ 100 Mq) EMP 
model can reproduce the observed colours of one line-offset 
LAE as well as the LAEs without offset (Figure 12a). A 
Pop III model (~ 1 Myr and ~ 100 Mq) can reproduce all 
but one LAEs discussed here (Figure 12b). The last LAE 
(object d) shows the largest offset between NB359 position 
and Lya, thus, the possibility of a foreground contamina- 
tion in NB359 is the highest. Finally, we conclude that the 
Pop III model is the most favourable stellar population for 
the LyC emitting LAEs although the EMP model may be 
still compatible. If we consider a combination of two stellar 
populations, 'primordial' and 'normal' sub-solar metallicity 
with dust, the UV two-colour diagram of all the LAEs and 
LBGs can be well explained (Figure 13). In this case, the 
mass fraction of the 'primordial' stars is estimated to be 
0.1-10% or more for LAEs with V - i' < 0.3, a few 0.1% 
for LAEs with V — i > 0.3 (but these may be foreground 
contamination), and 0.01-1 % for LBGs, depending on the 
age of the galaxies, dust content, and IGM attenuation (Fig- 
ure 13). We also note that the LBGs can be explained by 
normal stellar populations and EMP or Pop III stars are not 
mandatory for them. 

The Lyman 'bump' galaxies presented in this paper is a 
new population of high-z galaxies. Indeed, they were missed 
in standard 'drop-out' surveys because they are not 'drop- 
out' but have an excess flux in LyC: Lyman 'bump'. Because 
we started from a sample selected based on a narrowband 
excess by the Ly a emission line, we did find another e xcess 
in LyC. However. [Steidel. Pettini. fc Adelberger! (|200lh and 
IShaplev et al.l (|2006l ) started from a LBG sample, i.e. 'drop- 
out' sample, and thus, they could not find any galaxies with 
a Lyman 'bump'. This new population of galaxies probably 
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Phase-II: LAEs with LyC? 
(initial metal enrichment) 




Figure 14. Schematic picture of an evolutionary scenario of Ly- 
man 'bump' galaxies. 

often emits a strong Lya emission because they should in- 
trinsically emit strong LyC although a half of the LyC may 
directly escape. Therefore, starting from a LAE sample se- 
lected by a narrowband would be the best to sample the Ly- 
man 'bump' galaxies. On the other hand, if all the LAEs had 
Lyman 'bump', i.e. no 'drop-out', LAEs and LBGs would 
separate well. However, a si gnificant fraction of LBGs over- 
lap with LAEs in fact (e.g., IShaplev et al.|[2003l ; iNoll et all 
l2004l ; IVerhamme et al1l2008r ). This may imply that the frac- 
tion of LAEs having Lyman 'bump' is small although the 
IGM attenuation may play a role to hide the Lyman 'bump' 
and to make 'drop-out'. 

A caveat is the lack of the confirmation of the Lyman 
limit 'bump' in spectroscopy. Our argument of finding the 
Lyman limit 'bump' is only based on narrowband photome- 
try. We should confirm the feature in spectroscopy in future. 

6.1 A possible evolutionary scenario of Lyman 
'bump' galaxies 

The model with two stellar populations shown in Figure 13 
simultaneously reproduces the strength of the Lyman limit 
'bump' and non-ionizing UV slope of the LAEs and LBGs 
detected in LyC by 109. In the model, we have assumed 
that a primordial stellar population which is dust-free coex- 
ists with 'normal' stars with sub-solar metallicity and dust. 
We here discuss a possible scenario which may justify this 
assumption. Figure 14 shows a schematic picture of the sce- 
nario. 

The hierarchical structure formation scenario ar- 
gues that a galaxy is assembled from many subgalactic 



minihaloes. Cosmological hydrodynamics simulations have 
shown that the galaxy assembly occurs through cosmic fil- 
aments and galaxie s reside at nodes of the filaments (e.g., 
iBromm et alTl2009l ). iGreif et al.1 (120081 ) examined the first 
galaxy assembly at z ~ 10 and found that in a 5 x 10 7 
Mq halo, ~ 10 very massive metal-free stars, so-called Pop- 
ulation III. 1, were formed before the assembly. As a ra- 
diative feedback by the Pop III. 1 stars, pristine gas which 
was ionized once by the Pop III. 1 stars can be more abun- 
dant in H2 and HD than before, cool more efficiently, and 
form intermed iate-mass metal-free stars, so-cal led Popula- 
tion II I.2 (e.g.. lYoshida et al.ll2007l ) . As a result. IGreif et all 
(2008) expected that gas infalling onto the central part of 
a galaxy in the potential well likely forms Pop III. 2 during 
the assembly if the gas remains primordial. The Pop III. 1 
and III. 2 stars emit intense LyC which produces cosmolog- 
ical H 11 and He ill regions. The esc ape fraction of LyC is 
expected to be ver y high (60-90%; iKitavama et al.l |2004| ; 
lYoshida et aT1l2007D and the nebular LyC can make a sig- 
nificant Lyman limit 'bump' (Phase-I). 

While very massive Pop III. 1 may become blackhole di- 
rectly and make no contribution to the metal enrichment, 
Pop III. 2 stars will end their life with core-collapse super- 
novae (SNe) and enrich surrounding gas with metal (e.g., 
iHeeer et alj|2003l ). At t he same time, dust grains are also 
formed by the SNe (e.g-. lNozawa et afll200&t ). The metal and 
dust enrichment starts from the central part of the galaxy 
where the mode of star formation ch anges into the 'nor- 
mal' mode (e.g.. lSchneider et al.|[2006l ). If the multiple SNe 
blast waves make a supershell which contains dust grains, 
the dust attenuation becomes a screen geometry in which 
the dust attenuation law is determined only by dust optical 
properties. If the dust properties resemble the SMC which 
is also making an intense star formation, the reddening law 
has a steep wavelength dependence like the SMC extinction 
law which is assumed in the upper panel of Figure 13. On 
the other hand, primordial gas and Pop III stars may be 
still infalling onto the galaxy from the IGM. These infalling 
primordial stars can make a significant Lyman limit 'bump' 
(Phase-II). 

As the star formation proceeds in the main part of the 
galaxy, dust and metal are mixed well in the ISM. In such 
a medium, the wavelength dependence of the dust attenua- 
tion law beco mes weak and expected to be s i milar to the 
Calze tti law l|Fischera. Dopita. fc Sutherland! 120031 : llnoud 
120051 ) which is assumed in the lower panel of Figure 13. 
Even in this late stage, a small portion of the surrounding 
IGM may remain primordial if the metal enrichment in the 
IGM is not efficient. The primordial gas may form Pop III 
stars during its infalling process through a cosmic filament 
(Phase-Ill). 

With this scenario, we can explain the two-colour dia- 
gram as follows; the LAEs correspond to the Phase-II (or 
possibly Phase-I). If the age of the galaxies is 10-100 Myr, 
the primordial mass fraction is ~ 1-10% or more in some 
cases. The amount of the dust attenuation for the underly- 
ing population is E(B — V) = 0-0.3. The LBGs correspond 
to the Phase-Ill because of their relatively large stellar mass 
(~ 10 10 " 11 M e ). If the age is 100-1,000 Myr, the primordial 
mass fraction is < 1% and typically ~ 0.1%. The amount of 
the dust attenuation is E(B — V)= 0-0.3 which is similar to 
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LAEs but the attenuation laws may be different from each 
other. 

We can also explain another observational fact: spa- 
tial displacements of LyC and non-ionizing UV of the LBGs 
found in 109. As shown in Fig. 2 of 109, LyC detected po- 
sition coincides with a substructure found in the periph- 
ery of the main component of the LBGs, although some of 
these LyC emitting spots m ight be foreground contamina- 
tions (|Vanzella et al.ll2010al ). In our scenario, we expect that 
primordial stars emitting LyC exist in the periphery of the 
galaxies (see also iTornatore et al.l l2007). This is exactly the 
finding from the LBGs by 109. We also expect that LAEs 
show a smaller offset between LyC and non-ionizing UV than 
LBGs. According to 109, the offsets of the LAEs are not sig- 
nificant, and thus, smaller than those found in the LBGs. 
Again, this is consistent with our scenario. 

On the other hand, we have found another spatial dis- 
placement between continuum including LyC and Lya emis- 
sion line in some LAEs. While these continua may not come 
from the LAEs but from foreground objects, the probability 
that all the LAEs with line offset have a foreground contam- 
ination is very low as seen in §4.2.2. We have also proposed 
a possible mechanism to produce such a line offset in §4.2.3: 
the LAEs are composed of several unresolved substructures 
which have different Lya EWs. This mechanism may fit in 
our scenario as follows: in the Phase-II, the continuum in- 
cluding LyC comes from substructures in the periphery and 
Lya emission comes from the dusty central component, be- 
cause the Lya transmission can be higher than the contin- 
uum trans mission if the dust is confined in clumps in the 
supershell (|Neufeldlll99lh . 



6.2 Primordial stars at z ~ 3 and intergalactic 
metal enrichment 

Massive (~ 100 M Q ) EMP or Pop III stars exist in z ~ 3 
LAEs if our interpretation is correct. Since these primor- 
dial stars are short-lived (a few Myr) because of their mass, 
pristine gas should remain to form them until z ~ 3. 
This argument is not consistent with a scenario of global 
metal enrichment by Pop III stars at very early epoch 
(Mac kev. Bromm. fc HernquistJ 120031 ). in which the entire 
universe was enriched to a metallicity floor as Z > 1O _3 Z0 
by Pop III stars at z > 20. On the other hand, recent studies 
suggest inhomogeneous metal enrichment in the IGM and an 
extended epoch of Pop III star formation as described below. 

Observations of metal absorbers in high-z QSO spectra 
show that metal enrichment in Lya forest is inhomogeneous 
dSimcoe. Sargent, fc Rau ch 2004; Beck er. Rauch. fe Sargentl 
120091 ). For example. Simcoe et al. (2004) show that ~ 30% 
of lines of sight at z ~ 2.5 have metallicity lower than the 
detection limit of Z ~ W~ s Zq. Cosmological hydrodynam- 
ics simulations also show that metal enrichment in the IGM 
proceeds in a very inhomogeneous way and that a significant 
volume in the IGM rem ai ns pristine until a lower redshift 
1 Tornatore et al. 20071: Oppenheimer. Dave. fc Finlator] 



pristine mater ials kept until z ~ 3 (jTornatore et al.ll2007l : 
I Johnson! |20ich . In this sense, our argument of the discovery 
of EMP or metal-free stars at z ~ 3 is consistent with 
recent studies of IGM metal enrichment. 

The LAEs which may contain primordia l stars reside 
in a massive proto-cluster re gion, SSA22 field |Steidel et al.l 
1 19981: lHavashino et al.ll2004l ). This field is very peculiar in 
an enhanced activity of LAE s/LABs (M04 ), dusty sub- 
mm g alaxies (|Tamura et al.ll2~009l ) , and AGNs ( Lehmer et alj 
120091 ). Such activities may proceed IGM metal enrich- 
ment efficiently by metal-rich outflows wh i ch are observed 
ubiquitously around LBGs ( Ste idel et alj I2010T ). On the 
other hand, the LAEs ' clustering in this field is very low 
l|Havashino et all I2004T ) and this indicates that the LAEs 
do not lie in the high density peak of the proto-cluster 
but in the less dense periphery of the large-scale filaments 
IjShimizu et al.l 120071 ). The metal enrichment may not pro- 
ceed well yet there. However, the exact three-dimensional 
distribution of the LAEs and other types of galaxies by 
spectroscopy is not obtained yet and is an interesting fu- 
ture work. 



6.3 Implications for ionizing background and 
reionization 

The existence of massive EMP or Pop III stars in LAEs 
may have an impact on the studies of ionizing background 
and reionization. Based on the observed cosmic ev o lution 
of ionizing background, llnoue. Iwata. fc Deharvend (|2006T ) 
have suggested higher LyC emissivity (or escape fraction) 
of galaxies at higher-z. Such an evolving LyC emissivity 
is supported by galact i c LyC observation s |Shaplev et al.l 
l200fj; ISiana et al.ll2010l: iBridge et ai]|2010l ) and simulations 
iRazoumov fc Sommer-La rscn 200 dl2010h and is favourable 



20091: iTrenti. Stiavelli. fc Shulll 120091 : IWiersma et alj 1200a 
Maio et alj|2010l ). For example, Oppenheimer et al. (2009) 
show that only ~ 1% volume in the IGM is enriched to 
Z > W~' A Z Q by z = 5. Therefore, IGM metal enrichment 
is likely to proceed in a patchy way and we can expect a 
significant chance to form EMP or metal-free stars from 



for the cosmic reionization ( Bolton fc Haeh nclt 2007). On 
the other hand, the physical mechanism of the evolving LyC 
emissivity is still uncertain. Some radiative transfer sim- 
ulations suggest that l ess massive galaxies have larger es- 
cape fraction of LyC ([R azoumov & Sommer-Larscn |2010| : 
lYaiima. Choi, fc Nagamind |2010| ). IVanzella et alj (|2010d ) 
have found a possible luminosity dependence of the LyC 
emissivity from a large LBG sample. If such less massive 
galaxies are common at higher-z, the global LyC emissivity 
becomes larger at higher-z. Another possibility is that galax- 
ies at higher-z have a higher mass fraction of massive EMP 
or Pop III stars and their LyC emissivity is intrinsically 
higher. In §5.5, we have found that more than an order of 
magnitude larger primordial mass fraction in the LAEs (1- 
10%) than the LBGs (0.1%). If this larger primordial mass 
fraction is true for general LAEs and th e LAE population i s 
more common at higher-z as found bv lOuchi et al.l |2008), 
the galactic LyC emissivity becomes higher at higher-z. Of 
course, these two scenarios can be combined: the less massive 
LAE population has larger intrinsic LyC emissivity due to 
primordial stars and larger escape fraction. Indeed, to make 
a strong Lyman 'bump', we need a large escape fraction of 
~ 0.5 (Paper I). Therefore, such less massive Lyman bump 
galaxies may play a significant role in the cosmic reionization 
at z > 6. 
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